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A recent proposal seeks to explain the anomaly in the HERA high-Q^ 
neutral current e^p data by modifying the parton distribution functions at 
moderate-to- high-x and large . We investigate the consequences of this pro- 
posal for the HERA high-Q^ data itself, in the neutral current e~p channel 
and especially in the charged current e^p channels. We find that there are 
striking signatures in the charged current which already serve to rule out many 
possibilities, including the (symmetric) intrinsic charm hypothesis. For those 
not ruled out, interesting signals are predicted in e~p scattering. 
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1 Introduction 



Recently, the HI [||] and ZEUS collaborations at HERA announced an anomaly 
at high-Q^ in the e^p eX neutral current (NC) channel. Using a combined ac- 
cumulated luminosity of 34.3 pb~^ in e'^p — > eX mode at ^/s = 300 GeV, the two 
experiments have observed 24 events with > 15000 GeV^ against a Standard 
Model (SM) expectation of 13.4 ± 1.0, and 6 events with > 25000 GeV^ against 
an expectation of only 1.52 ± 0.18. Furthermore, the high-Q^ events are clustered at 
Bj or ken- a; values near 0.4 to 0.5. 

A number of authors have presented proposals for new physics which might explain 
the NC anomaly, including: new contact interactions |, s-channel leptoquark 
production S ^, and R-parity violating supersymmetry 0, |^, and related propos- 
als P). There exist a variety of criticisms for each of the proposed ideas, from their 
highly speculative nature, to very concrete flavor problems which they all share 0. 
As such, another somewhat less speculative idea has been proposed in which the 
anomaly in the HERA data is simply further evidence of our inability to calculate, 
and often even reliably fit, the parton probability distribution functions (PDFs) re- 
sulting from the non-perturbative dynamics inside the proton. That is, Kuhlmann, 
et a/ [0, have suggested a way by which the behavior of the PDFs at large-x can be 
modified to explain part, or all, of the HERA NC data without disrupting the fits of 
the old low-Q^, low-to-moderate-x data. 

It is actually a simple exercise to show that increasing the parton densities at 

^ 20, 000 GeV^ and X ~ 0.4 to 0.5 can in fact fit the NC data. It is a far more 
complicated question whether these changes are consistent with all other world data. 
We will not consider here the validity of the claim that such a fit can be done which 
is consistent with the low-Q^ data. Instead we will show that there can be dramatic 
consequences for the HERA data itself in the charged current mode at high-Q^. We 
will find that the HERA data already rules out a number of theoretically acceptable 
scenarios for modifying the PDFs to agree with the NC data. 

HERA is capable of running in two modes: e~p and e~^p. In the former mode 
HI and ZEUS have accumulated 1.53 pb~^ of data but have observed no statistically 
significant deviations from the SM. Further, the experiments differentiate between 
final state eX and uX, where the neutrino is identified through its missing pt- HI 
has also announced its findings in the e~^p vX charged current (CC) channel. 
They find 3 events at > 20000 GeV^ with an expectation of 0.74 ± 0.39, but no 
events with Q'^ > 25000 GeV^. ZEUS has not announced its CC data as of this date. 
Although compared to the NC data, the present CC data is much sparser, one can 
still conclude from it that there can be no deviations from the SM by more than 
a factor of 2 or 3 in that channel. Whether there is in fact any deviation at all is 
still too uncertain to say. (Explanations of the possible CC excess involving non-SM 
physics have been considered by [|10]-) However, given the sizes of the effects which 
we are going to find, more data in e^p, and especially e~p, mode will provide very 
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strong constraints on the viability of this suggestion. 

2 The proposal 

The class of proposals considered by Kuhlmann, a/ 0], for explaining the HERA 
anomaly can be thought of as variations on the so-called "intrinsic charm" scenario |]TT| 
in which it is posited that there is some non-perturbative, valence contribution to the 
charm quark distributions at — > 0. (Current fits assume that there is no valence 
charm component in the proton and generate non-zero contributions at higher 
only through perturbative renormalization group flow, i.e., gluon splitting. In this 
scenario one can think of the valence structure of the proton being uudcc versus the 
usual uud.) 

Though the underlying dynamics (and motivation) may be subtle, the proposal 
itself is straightforward: increase by hand the parton densities inside the proton 
at X ~ 0.5 and large to the point that the NC cross-section of the Standard 
Model matches that observed at HERA. Such an effect can arise naturally from non- 
perturbative dynamics at low if the dynamics produce a narrow "bump" in the 
PDFs at low Q2 and very large a; ~ 1. This bump would migrate down to lower x 
as one flows (through the renormalization group) up to higher Q^. Such a structure 
is difficult to rule out; the data at low and very high x is limited and extractions 
of the structure functions are problematic due to non-perturbative and higher-twist 
effects which can be important at large x. Fits available now typically use only data 
with X < 0.8. The proposal in emphasized enhancements to the w-quark density 
rather than that of the c-quark, but the difference is irrelevant from the point of view 
of this paper because electroweak physics does not distinguish among the generations. 
Enhancement of any or all of the parton densities could in principle explain the HERA 
data. 

This proposal has strong advantages and disadvantages. In its favor, it requires 
no new physics beyond the Standard Model and thus automatically solves the flavor 
problem associated with many of the new physics interpretations ||^. To its detri- 
ment, it invokes non-perturbative QCD effects which are not calculable and cannot 
be predicted ah initio, and its consistency with the moderate data taken at both 
HERA and the Tevatron has not been fully studied. 

3 High-Q2 tests at HERA 

Putting aside any questions of how well such a proposal can really do at explaining 
the HERA NC data while remaining consistent with all other world data, one can 
ask: what are the consequences of such a suggestion for the NC process in e~ mode, 
and for CC processes in either mode? 

The NC result is simple: in this scenario, the NC cross-sections in both e~p and 
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e+p modes scale by approximately the same amount. The Z couplings introduce a 
small helicity dependence that keeps the two modes from being exactly the same. 
However, the differences between the two are less than 10% because the photon con- 
tribution typically dominates the total cross- sect ion. Thus this proposal predicts 
that once HERA has accumulated enough data in e~p mode, they will observe a NC 
anomaly there as well, and it should be roughly the same size as that observed in 

In discussing the CC predictions, a number of complications arise, most stemming 
from the fact that the x and dependences of the NC and CC cross-sections are 
somewhat different. Nonetheless, we will present fairly precise heuristic arguments 
about the sizes of CC effects, which we will then check in a full numerical calculation. 

One key simplification is this: in the SM, even at high-Q^, the NC scattering is 
largely dominated by virtual photon exchange in the t-channel. This need not have 
been so, since at >> vtiz there is no additional kinematic suppression of the Z 
contributions; however, the Z coupling to quarks is generally smaller than that of 
the photon. For the arguments that follow, we will ignore the Z contributions and 
reintroduce them only when going to the full calculation in the next section. 

For photon exchange alone, the NC cross-section at high-Q^ behaves as: 



1 dOT^c , . , 
; OC U\X) + 

X ax 



u{x) + J [d{x) + d{x)} + ■ ■ 



where u{x) is the w-quark parton probability distribution function (PDF) inside a 
proton, u[x) is the w-antiquark PDF, etc., and the ellipses represent heavier (s, c, b, 
t) quarks. The factor of 1/4 is the relative charge-squared of u- and c?-type quarks. 

There are four orthogonal classes of changes to the PDF's that can be considered, 
each corresponding to enhancing the densities of either u, u, d or d individually. That 
we do not have to consider changes to the charm and strange densities is clear, since 
the physics in question cannot distinguish u from c, or d from s. (This is not a general 
statement for all experiments at all Q^. For example, it does not hold even at HERA 
if HI and/or ZEUS could tag prompt charm production.) Each of these four cases 
leads to a distinctive CC signature. 

We will parameterize the effect of an "intrinsic quark" component on the PDFs 

by: 

= Qo{x) + qintix) = qoix)eg{x) (2) 
where q{x) is the total parton distribution, qo{x) is the usual fit distribution (we 



will use the CTEQ3 set whenever we have to make an explicit choice [|12[, similar 



results would follow from the MRS sets [0), qmt{x) is the intrinsic component and 
^q{^) ^ 1 parameterizes the effects of the enhanced component. (The dependence 
of q{x) and e{x) is implicit and of little relevance to what follows since we will only 
consider scattering within a small range of values; therefore the renormalization 
group running with Q"^ can be ignored.) For simplicity, suppose that eq{x) is exactly 
1 everywhere except in a small range of x centered on x = xq for which it is much 
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greater and roughly constant {i.e., it is roughly a top hat distribution): 



eg{x) = 1 + r]g9{x - xo + S)9{xo + 5 - x) (3) 

where r]q > 0, 6{x) is the usual step function, and 6 is the half-width of the top 
hat. Though the numerical arguments do not rely on making this assumption, it does 
make the pedagogy simpler. 

To begin, suppose that the correct fit to the HERA NC data is achieved by 
changing only u{x) such that Xq, 6 and e^lxo) have certain fit values. This will be 
the canonical scenario to which wc compare all others. For example, suppose that 
instead of changing the u{x) PDFs, we would like to do a fit for which only d{x) is 
changed. Then to produce the same NC cross-section that e„(xo) provided in the 
It-case, ed{xo) must shift by a larger amount, though at the same x — xq. That 
the shift must be larger is clear, because rf-quark effects on the NC cross-section are 
suppressed both by electric charge (the 1/4) and by the smaller li-quark content in 
the proton, do{xo)/uo{xo). In terms of an equation, setting the shift due to the in 
the one case equal to the new shift by e^: 

1 1 

eu{xo)uo{xo) + -do{xo) = uo{xo) + -ed{xo)do{xo), (4) 

giving = 1 + 4(eu — l)(-Uo/c?o), with all functions evaluated at x = Xq. To get the 
same effect by enhancing u{x) alone requires eu = ! + — 1)(mo/mo) a.t x = xq. 
Similarly, changing d{x) alone is identical to the case of d{x), but with d — > d in the 
expression above. 

What effects do these changes have on the charged current? In terms of the PDFs, 
the CC cross-sections scale as: 

1 rlrr-^ 

(X u{x) + d{x){l - yf + ■ ■ ■ (6) 



X dxdy 

where the +[— ] superscript denotes scattering in e'^p[e~p\ mode, and y has its usual 
definition in deep inelastic scattering. In the center-of-mass frame, y — sin^(i?/2) 
where i) is the scattering angle; its presence in the expressions above follows 
trivially from angular momentum conservation in the {V — A) scattering process. 
Note also that for the usual PDFs at large x, to a very good approximation a^c ^ 
d{x){l — yY and (t^c u{x). 

For the moment, in order to examine the effects of enhancing the PDFs on the 
CC differential cross-section, we will restrict ourselves to the region around 
of half-width 5. We will denote the differential cross-section in this small region 
daQQixo) as a shorthand. In the first case discussed above, in which only u{xq) is 
changed in response to the NC data, the CC signal in e+ mode is unchanged, while 
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daQ^{xo) increases by eu{xo)- Thus the same relative change in the NC data will also 
occur in the e~p CC data, at least in the neighborhood of x = Xq. However, because 
the kinematic dependence on x (at high-Q^) is roughly the same in e~^u e^u as in 
e~M ud, the relative scaling of the NC to CC in the xo-region will hold for all x. 
(The preceding statement is exact in the limit ^ and when the Z contribution 
to NC can be ignored. Since these two conditions are not simultaneously satisfied in 
general, there will be important corrections to the results of these heuristic arguments. 
Nonetheless, the results derived here do not differ greatly from the exact results, as 
we will see.) Thus, to have a concrete example, if e„(xo) is chosen to double the high 

NC cross-section, it will also (approximately) double the high-Q^ CC cross-section 
in e~ mode, while having no effect in e"*" mode. 

The effect on the CC is more marked in the case of using d{x) to explain the NC 
anomaly. As we said above, to have the same effect on the NC data as e„, the must 
be roughly 4:[uo{xo) / do{xo)] times bigger than the corresponding e„ would have been. 
And since to a good approximation 

d^cci^o) oc ed{xo) ~ 4 ^°|^°| e^(xo), (7) 

do{xo) 

it will scale by the same amount. Consider again our example, where we demanded 
that e„(xo) double the high-Q^ NC data. For xq ~ 0.5 one finds uo{xo) ~ ido^xo). 
Then to fit the same NC data, 6^(2:0) must be 4 ■ 4 = 16 times larger than eu{xo). 
Moving this enhancement to the CC one finds that da^^^lxo) increases by a factor 
of 16e„(a;o), while a^Q is unchanged. Again the explicit kinematic dependences on x 
are (approximately) the same in the CC and NC, so the overall a^f-; will increase by 
roughly a factor of 32 in the high-Q^ region. 

The same arguments go through for u and d. For u, da^^^Xo) scales by a factor 

d.,,(x.)ocl+|^-=.|^-, (8) 



while is unchanged. For explaining the HERA anomaly we would be interested 
in ~ 2 X 10'' GeV^ and x ~ 0.5, and thus y = Q'^/{xs) ~ 0.5. Then in our 
recurring example wherein the NC signal is doubled, daQf-;{xo) scales by a factor of 
I6eu, leading to an overall scaling of a^Q again by 32. 

In the final case, with enhanced d{x), daQQ{xo) scales by a factor 

d(T^c(xo) ocl + 4(e,-l)(l-y)2 (9) 



while aQ<^ is unchanged. In our recurring example, with y ~ 0.5, we find da^^lxo) 
scales by approximately e„, so that the full a^^ at high-Q^ is expected to be double 
the usual SM prediction. 

The various scalings of the CC mode are summarized in Table ^ For each parton, 
if we assume that it alone must have its density scaled to explain the HERA NC data. 
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Part on 


NC Factor 


CC+ Factor 


CC" Factor 


u 

CL 

u 
d 


i + 4(e„ - ij^ 
l + 4(e„-l)| 


1 

^ J- Aft: 1 M 

^ {1-yy do 
1 


1 

-L 

1 

l + 4(e„-l)(l-i/)2 



Table 1: For each parton is shown the relative sizes of contributions needed to explain 
the HERA NC data, as well as the predictions for the enhancement of the CC signal, 
daQ(j{xo)- All quantities are implicitly functions of x and are to be evaluated at 
X = Xq. These results are derived in the heuristic scenario discussed in the text. Note 
that for Ca ^ 1 the scaling factors are independent of e„. 

the size of the enhancement needed (relative to the enhancement e„ for the w-quark) 
is given in the second column. In the third and fourth column are then shown the 
relative enhancements of daQ(j{xo) given a fit to the NC. 

Some profound results can already be extracted from these simple considerations. 
If the PDF's of (m, d, u, d) are each changed respectively so as to double the high-Q^ 
NC data, then the e"*" CC data will increase by factors of about (1, 32, 32, 1) while 
the e~ CC data will increase by (2, 1, 1, 2). Since doubling the NC cross-section at 
high-Q^ is roughly consistent with the HERA data, we can conclude that we either 
expect extremely large enhancements in the e"^ CC signal, or none at all. The data 
from HI is not consistent with an extremely large enhancement (such as a factor of 
32), so we can conclude that: (1) changing the PDF's of a ci-type quark or a w-type 
quark to explain the HERA NC data is inconsistent with the HERA CC data; (2) 
any other explanation involving modifications to the w-type or c/-type PDFs will not 
lead to a CC signal in mode, but will have one in e~ mode; (3) if HERA sees 
an enhancement in the e^p CC channel comparable to the one in the NC channel, 
more than one PDF must be modified. Note that (2) above does not preclude a small 
CC enhancement in mode; but any such enhancement will not be enough to also 
explain the NC data. 

Up until now, our arguments have ignored the complications introduced by keeping 
the Z contributions in the NC, and by the non-zero W mass in the CC. To show how 
these results are modified in a complete calculation, we have numerically evaluated 
the CC signals which would be induced by fitting to the NC anomaly. In Figure |l|, 
we have considered each of the u, d, u and d cases separately as a solution to the 
HERA NC anomaly. (We are only using the HI NC and CC data since as of this 
writing ZEUS has not announced their CC results.) We have chosen in each case 
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1 2 3 4 5 



12 3 4 



Q2 (104 QeV2) 

Figure 1: Differential cross-section for CC e^p and e~p scattering as a function of 
Q^. Dotted lines are the SM prediction with unmodified PDFs; solid lines are with 
modified PDFs. 

for eq{x) to have a top-hat form as in Eq. (j^) with xq = 0.45 and S = 0.05. This 
region in x envelopes the bulk of the HI high events. The overall normalization 
rjq is chosen to provide the best possible fit to the HI data. (The best fit values 
of rjg for q = u,d,u,d were found to be 7.5, 95, 640, 450 using the CTEQ3M PDFs. 
rjg for the sea quarks shows a strong dependence on the PDFs used, because of the 
uncertainty in the sea quark distributions at moderate x and large Q^. However the 
overall enhancement of the CC signal is independent of the choice.) Then we have 
plotted the resulting da'^fj/dQ'^ for each possibility. 

The Figure is divided into two frames for the CC processes e^p VX and 
e~p vX separately. For the e^p mode, the enhancement of the CC cross-sections 
is marked, consistent with our heuristic derivation of an enhancement factor of ~ 32. 
For the e~p mode, the enhancements are much smaller, once again consistent with 
our expectation of factors of ~ 2 only. 

Notice from the figure that the peak enhancements as a function of can be 
much larger than those of the total integrated cross-sections, which we give in Table ^ 
There we have done the integration in the Figure, for > 10, 000 GeV^ and > 
20, 000 GeV . The most important HI cuts have been included in the calculation: 
y < 0.9 and PT,miss > 50 GeV. Our pedagogical derivation of the enhancements has 
been shown to work reasonably well, though not exactly. The two PDFs which affect 
the e~^p mode both induce very large corrections in the CC if they are invoked to 
explain the NC data. However the two PDFs which affect the e~p mode induce 
small, but observable, corrections. Given high statistics, HERA should be able to 
probe even the hardest case, that of changing d{x) to explain the NC data. With 
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PDF 


Mode 


> 10, 000 GeV 


> 20, 000 GeV 


u 


e~p 


1.9 


2.9 


d 


e^p 


10 


21 


u 


e+j9 


5.5 


25 


d 


e"p 


1.3 


1.2 



Table 2: Multiplicative enhancements of the CC cross-sections for above the 
indicated value, in the mode indicated. The PDFs in the first column have been 
changed to give a best fit to the HI NC data. 

the current data, it already appears that altering d{x) or u{x) cannot be invoked to 
explain the NC anomaly. 

The possibility exists of probing more complicated combinations of modifications 
because the cross-sections at HERA are linear in the PDFs. Therefore in a scenario 
in which several of the PDFs are modified so that the total modification is a sum of 
individual parton modifications, weighted by some then the resulting CC signals 
are also linear combinations of those shown, again weighted by the a^. For example, 
if enhancements to the u and d are such that each one contributes half of the NC 
excess, then the CC signal will be an average of the individual signals for u and d. In 
particular, the intrinsic charm scenario, for which ec = e-c is responsible for explaining 
the HERA data, the da(^(j{xQ) is scaled by 1 -|- (e^ — l)(Mo/2(io)/(l — ?/)^ and da'^fji^xo) 
is scaled by (l-t-eu)/2 where = 1 -|-2(ec — 1)(co/mo)- This linear behavior provides a 
powerful, and simple, way to disentangle the form of the modifications given the size 
of the excesses over SM in the two CC modes. Unfortunately, with only 2 data points 
for 4 unknowns, one cannot deduce the full answer using CC data alone. However, 
it is already enough to rule out the intrinsic charm scenario as an explanation of the 
NC data, given the numbers in Table ^ 

4 Conclusions 

Modifying the PDFs to explain the HERA e^p NC data implies an immediate modifi- 
cation (by about the same relative size) of the e~p NC data, and also implies striking 
patterns of modification in the CC data. Such modifications should be easy to ob- 
serve given the current size of the NC anomaly. Further, they already rule out the 
intrinsic charm scenario (with equal modifications to c(x) and c(x)), or any other 
attempt to use modifications of d-type or w-type distribution functions to account for 
the anomaly in the HERA NC data. Our results for the CC channel are summarized 
in Table ^ and Figure |l|. 

If an excess in the CC signal is detected at HERA in the e^p data of the size 
currently suggested by HI, then one must go to a scenario in which more than one 
PDF are modified, but such that one dominates the contribution to the e^p NC 
excess. This will lead to an interesting signal at HERA in the e^p CC mode. In 
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particular, once significant luminosity has been collected in e~p mode at HERA, one 
should also be able to probe solutions to the NC anomaly which involve only changing 
the distributions of u- and rf-type quarks. 
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